Dark Matter distribution in the Milky Way: microlensing and dynamical constraints 
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We show that current microlensing and dynamical observations of the Galaxy permit to set 
interesting constraints on the Dark Matter local density and profile slope towards the galactic centre. 
Assuming state-of-the-art models for the distribution of baryons in the Galaxy, we find that the most 
commonly discussed Dark Matter profiles [viz. Navarro-Frenk- White and Einasto) are consistent 
with microlensing and dynamical observations, while extreme adiabatically compressed profiles are 
robustly ruled out. When a baryonic model that also includes a description of the gas is adopted, 
our analysis provides a determination of the local Dark Matter density, po — 0.20 — 0.56 GeV/cm 
at Icr, that is found to be compatible with estimates in the literature based on different techniques. 



I. INTRODUCTION 

The presence of large amounts of Dark Matter (DM) in 
the Universe, and in particular in the Milky Way, is estab- 
lished on sound observational grounds [TH3]. However, a 
detailed description of the Dark Matter distribution in 
the Galaxy is hard to achieve, despite the tremendous 
progress in numerical simulations of galaxy-sized objects 
over the last few years m?] , therefore precluding, among 
other things, a precise interpretation of direct and indi- 
rect Dark Matter searches (see e.g. Refs. [8Hl3]). In order 
to constrain Milky Way mass models, several observables 
have been used in the literature, including star counts, 
the motion of gas and stars or microlensing events. Here, 
we focus on two key probes, viz. microlensing observa- 
tions and dynamical measurements, and show that in- 
teresting constraints can be set on the DM distribution 
using those measurements only. This provides comple- 
mentary evidence for the existence of Dark Matter in the 
Galaxy, with a distribution compatible with that inferred 
from other observables. 

Gravitational microlensing has long been adopted as 
a tool to study the structure of the Galaxy [Ml [15]. 
This technique was suggested in [14] to probe the ex- 
istence of Massive Compact Halo Objects (MACHOs) by 
using source stars in the Magellanic Clouds, M31 and 
M33. Soon after the observation of the first microlens- 
ing events in the early 90s [TMT5] . it became apparent 
that MACHOs could not be a dominant component of 
the galactic dark halo [T^ - therefore a type of dark, 
non-compact matter was needed. Meanwhile, microlens- 
ing of stars towards the galactic centre proved useful in 
tuning the existent bulge, bar and disk models. Until to- 
day, several thousands of microlensing events have been 
collected - most notably by OGLE [HI [201 El] , MACHO 
[ia[21[23], EROS [nUMMS] and MOA [17] campaigns 
(see [2H1 ES] for a review) - allowing for precise estimates 
of the microlensing optical depth towards the galactic 
bulge, spiral arms and Magellanic Clouds. 

Another well-known tool to constrain the different 
components of our Galaxy is the rotation curve. Un- 



like microlensing - which is sensitive to the distribution 
of compact matter along the line of sight -, dynamical 
measurements constrain the total mass distribution. Sev- 
eral decades ago, the observation of flat rotation curves 
in external spiral galaxies provided convincing evidence 
for Dark Matter, but actually the rotation curve of the 
Milky Way (a spiral itself) is far from being precisely de- 
termined given our peculiar position within the galactic 
disk. In fact, although hosts of data have been gathered 
over the years [30], allowing for ever more accurate dy- 
namical models, sizeable uncertainties remain on Milky 
Way mass models, a circumstance that in turn affects our 
capability to constrain the DM distribution, as we shall 
see below. 

In a pioneering work 31: . Kuijken pointed out that a 
non-zero microlensing optical depth towards the galactic 
bulge sets a lower limit on the enclosed baryonic mass 
and thus a lower limit on the corresponding circular ve- 
locity within the solar circle. Using this result, Binney 
& Evans j32] derived an upper limit on the dark mat- 
ter contribution to the rotation curve, which led them 
to the conclusion that "the cuspy haloes favoured by the 
Cold Dark Matter cosmology (and its variants) are incon- 
sistent with the observational data on the Galaxy" [32j . 
This claim was however based on a preliminary determi- 
nation of the microlensing optical depth by the MACHO 
collaboration [33] (later replaced by a less constraining 
measurement from the same collaboration [53]) as well 
as a very simplified treatment of the galactic rotation 
curve. 

Here, we revisit this claim in light of recent microlens- 
ing observations towards a variety of galactic regions and 
up-to-date measurements of our Galaxy's rotation curve. 
A very significant improvement of our work with respect 
to the existing literature is a proper inclusion of all ex- 
perimental uncertainties regarding both microlensing and 
dynamical data, which makes our conclusions sound on 
statistical grounds. Furthermore, wc shall also use a wide 
range of state-of-the-art models for the galactic bary- 
onic component, so that modelling uncertainties are ap- 
propriately bracketed. Our ultimate aim is to combine 
microlensing and dynamical observables to draw robust 



constraints on the Dark Matter distribution in the inner 
region of the Milky Way. 



relevant bulge model is the one introduced by Zhao in 



II. MILKY WAY MODELLING 

In order to extract useful constraints on the DM dis- 
tribution from microlensing and dynamical observables, 
we need to specify a mass model for our Galaxy. In this 
work, we shall implement a bulge/bar region, a disk and 
a Dark Matter halo. The Sun sits on the galactic disk at 
(x0, 2/0,^0) = (i?0:0,0), where Rq is the galactocentric 
distance^. We consider the range Rq = 8.0 ± 0.5 kpc 
throughout the analysis, taking 8.0 kpc as our fiducial 
value. This encompasses recent ^S] and less recent [35] 
determininations as well as the 1985 lAU recommended 
values [37] . All scale distances in this Section are rescaled 
as Ro/8 kpc. 

The bulge at the centre of the Milky Way has long 
been suspected to be bar-shaped. This was confirmed 
by near-infrared observations |38j and the distribution 
of red clump giant stars [SSI HU], which produce con- 
vincing evidence for a triaxial bulge with its near end at 
positive galactic longitudes £ > and its major axis in- 
clined a& ~ 25° with respect to the galactic centre line of 
sight. Several families of triaxial mass density distribu- 
tions have been tested against the data, and two param- 
eterisations appear to work particularly well [33 [35] : 



Pbix,y,z) = po^te ''' 
Pb{x,y,z) = po,be~^ 
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where (a;', y' , z') are the coordinates along the major, in- 
termediate and minor axes, respectively, and the param- 
eters of these fitting formulae are [SS] 



E2: ab = 23.8°; {xb, yb, Zb) =(0.899, 0.386, 0.250) 
G2: ab = 24.9°; (x^, yb, Zb) =(1.239, 0.609, 0.438) 



Rq 
8 kpc 

Rq 
8 kpc 



kpc 



kpc . 



For both the E2 and G2 models, we implement an ex- 
ponential cutoff proportional to exp(— (_R' — i?m)^/2r§) 



(with i?'2 = a;'2 + y'- 



ro 



0.5 kpc) beyond the co- 



rotation radius Rm = 3.5 kpc as in [3H1 EI]. Another 



^ Actually, the Sun lies somewhat off the central plane at zq < 
30 pc 1341 . but this is a very small displacement compared to 
other quantities at play and thus we shall take the reasonable 
assumption zq = 0. 



Pb{x,y,z) =po,b s, 



-1.85, 



(Zhao) , (3) 



with si = {ql{x^ + y"^) + z'^)/zl, which is particularly 
suited to produce a steep rise in the rotation curve at 
r < 0.5 kpc. Inspired by j43) we shall consider an ax- 
isymmetrised version of this profile with Xb — yb — 
0.9 kpc(i?o/8 kpc), Zb = 0.4 kpc(i?o/8 kpc) and qb = 0.6. 
Finally, the recent and sophisticated model of Gardner 
et al [44J 115 is also implemented - in this case there are 
separate bar and bulge components: 

Pb{x, y, z) = /o,& {pbar{x, y, z) + pbuige{x, y, z)) (4) 

(Gardner et al) , 

where the bar density is given by a Ferrers model [46] 

(^ ,, ^\-i Po.bar (1 - rl)'^ for n < 1 
r{x,y,z)-< ^ forri>l ' 



Pba 



with n = 2, ab = 25°, Xb — 3.5 kpc(i?o/8 kpc), yb = 
1.4 kpc(i?o/8 kpc), Zb = 1.0 kpc(i?o/8 kpc), Mbar = 
10^° M0 ^44], and the bulge component is given in terms 
of the gravitational potential (pbuige defined in [45] . The 
mass density associated to 4>buige is readily obtained us- 
ing the Poisson equation. The constant /q j, is just a nor- 
malisation and its role will become apparent in Section 

imi 

As for the stellar disk in our Galaxy, it is well-known 
that there are a thin and a thick components correspond- 
ing to distinct star populations. Instead of covering a 
whole range of different vertical and radial disk profiles, 
we adopt two extreme cases: a thin-only disk and a 
thin-|-thick disk. Although the former model is somewhat 
unrealistic, it will prove useful in exploring the present 
uncertainty associated to disk modelling. For an alterna- 
tive disk model see Ref. i47j . Following [HI US] , 



Pd[R,z) = e H 



(1 - /3)sech^ -— +I3e'^^ (5) 
■qhi 

(Han & Gould) 



with 



R 



'2 - a;2 . 



y\H = 2.75 



Rn 



8 kpc 



kpc. 



hi = 0.270;:^^ kpc, h2 = 0.440-^ kpc. 



8 kpc 



tjiR) — max 



0.670,0.114- 



kpc 
R 



9.025 kpc J 

thin: /? = 0, po,d =4.4 x 10^/(i?o) M0/kpc^ 
thin+thick: /3 = 0.565 , po.d = 4.93 x lO'^ M0/kpc^ . 

Again, we shall also refer to the model of Gardner et al 
where both stellar and gas disks are modelled through 



the corresponding gravitational potentials (see details in 

05]). 

In the remainder of the work, five models for the galac- 
tic baryonic component (i.e. bulge/bar and disk) will be 
used: 

• Model 1: E2 bulge and thin+thick disk; 

• Model 2: G2 bulge and thin+thick disk; 

• Model 3: G2 bulge and thin disk; 

• Model 4: Zhao bulge and thin disk; and 

• Model 5: Gardner et al bulge/bar, stellar disk and 
gas disk. 

We emphasise that these fiducial setups reasonably 
bracket the uncertainties in modelling the galactic bary- 
onic component, and thus will be useful in assessing how 
our results depend on such modelling. 

The last piece missing in our Milky Way mass model is 
the Dark Matter halo. In view of the findings of numeri- 
cal simulations [H [3] , we implement spherically symmet- 
ric generalised Navarro-Frenk- White (NFW) and Einasto 
profiles: 



PDM{r) 
PDM{r) 



p,(r/r,)-"(l + r/rs)-3+" (NFW) (6) 



PsGXp 



- I 



(Einasto), (7) 



where r^ — x'^ + y'^ + z"^ , rg is the scale radius, ps is the 
scale density and a is the inner slope for the NFW profile 
and a shape parameter for the Einasto profile. Typical 
ranges found in A'^-body simulations are 0.9 < a < 1.2 in 
the case of NFW [HHj and 0.12 < a < 0.22 in the case 
of Einasto |3nHS2- In the following, the normalisation 
of the DM profile will be set by the local DM density, 
Pq = Pdm{Ro)i that can be easily cast in terms of Ps- 
We shall take a fiducial interval p^ = 0.4 ± 0.1 GeV/cm 
in line with recent determinations [131 [55H55] . 

Regarding the scale radius r^, DM-only simulations 
tell us that the virial concentration Cmr = fmrlfs of 
an object of mass Myij. = 10^^/i~^ Mq lies in the range 



logi 



0.9 - 1.1 (at Icr) 



which translates into 



Ts ^ 24 — 38 kpc. However, works where galactic dynam- 
ical observables are used to constrain a Milky Way mass 
model [33] (see also [57]) seem to favour higher concen- 
trations, Cmr = 10 — 25, or Ts ^ 12 — 30 kpc. Hence, we 



-1-15 
-10 



shall take in the following a rather wide range r^ = 20 
kpc. 

We stress that the Dark Matter profiles found in nu- 
merical simulations deviate significantly from spherical 
symmetry: while DM-only simulations lead to very pro- 
late shapes, the inclusion of baryons leads to more oblate 
(but still triaxial) DM distributions 58, 59J. In order to 
check the relevance of a non-spherical DM halo for our 
purposes, we shall also consider the NFW profile in equa- 
tion ^ with r replaced by m = {x^ + y'^ + z^ jq^)^!'^ and 
axis ratio q — 0.7. This represents an oblate profile with 
a shape compatible with the results of [54) . 



As of today, it is not clear how baryons affect the DM 
profile in our Galaxy. One possibility is that when the 
Galaxy formed and the baryons contracted towards the 
centre, the dissipationless component was dragged lead- 
ing to a steepening of the DM profile. This was proposed 
long ago in Ref. [60], where a simple adiabatic contrac- 
tion model was introduced (more refined models have 
also been constructed, see Refs. ^TJIS^): given a spher- 
ically symmetric initial total mass distribution Mi(< r^) 
and a final baryon distribution M(,(< r), the final DM 
distribution M£)m(< t) obeys [30] 



r (Mt(< r) + Mum{< r)) = r,Af,(< r,) 



(8) 



with Mdm{< r) = (1 - fb)Mi{< ri), fb being the to- 
tal baryonic fraction in the Galaxy. Taking Mi{< ri) to 
be the initial mass distribution associated to the NFW 
profile (p| and Mi,{< r) the baryonic mass distribution 
given by models 1-5, we can easily solve ([8]) for r while 
fixing ri and /;,. This sets the final mass distribution 
Mdm{< t) and (upon derivation) the DM profile. For 
the sake of completeness, in our analysis we shall also 
use the described procedure to model the DM profile and 
therefore test adiabatic contraction models against galac- 
tic dynamical observables. 

Finally, note that in the following, given a baryonic 
model and assumed rs, we let vary po and a. These will 
be the main phenomenological parameters in our analy- 
sis. 



III. MICROLENSING 

In this Section we present a brief overview of the grav- 
itational microlensing formalism and of the most recent 
observations. 



A. Theoretical framework 

Gravitational microlensing is a direct consequence of 
Einstein's General Relativity: as a massive compact ob- 
ject - the lens - moves across the line of sight towards 
a given luminous source, the source light gets deflected 
and two distinct images form, assuming a Schwarzschild 
lens. In the microlensing regime these images are not 
resolved and hence the observable effect is a temporary 
magnification of the source. For a point-like source and a 
point-like lens, the magnification of the source reads [14] 



A{t) 



u{tf 



u{t)./^^ 



(9) 



where u{t) is the transverse distance between the lens and 
the line of sight measured in units of the Einstein radius 
Re = ^i4GMi/c'')Di{l-Di/D,). Here, Di (D,) is the 
distance between the observer and the lens (source) and 
Ml is the lens mass. The typical duration of a microlens- 
ing event is set by the Einstein radius and the lens ve- 
locity in the plane perpendicular to the line of sight vt'- 



tE — Re/vt- For galactic lens objects with masses rang- 
ing from 10~^ M0 to 10^ M0, the typical microlensing 
time scale varies from a few hours to a few years [H] . 
Therefore, a thorough monitoring of galactic stars in dif- 
ferent directions can be used to study a whole population 
of low-luminosity objects along the line of sight, which 
would otherwise remain undetected. 

The probability of observing the microlensing of a 
given luminous source is driven by the mass distribu- 
tion of lenses pi along the line of sight. A particularly 
interesting quantity is the optical depth r which quanti- 
fies the probability that a lens comes within one Einstein 
radius of the line of sight, or equivalently that the magni- 
fication of the source exceeds A = 1.34 (cf. equation ^ 
with u = 1). If all sources towards galactic coordinates 
{£, b) are placed at the same distance Dg (we will relax 
this simplifying assumption for a more general formula- 
tion later on), the theoretically expected optical depth 
can be written as: 



Ti£,b,D,) 



AttG 



Ds 



dDipi{e,b,Di)Di[l 






(10) 

where pi refers to all possible lens objects (in our case 
the bulge/bar and the stellar disk). For our purposes 
it is important to stress that (i) microlensing is caused 
only by massive compact bodies, and not by gas nor Dark 
Matter; and (ii) the optical depth t is independent of the 
lens mass function, since the surface of the Einstein disk 
is proportional to Mi. As a consequence, measurements 
of the optical depth depend only on the line of sight in- 
tegral of the mass density of lenses pi as emphasised in 
equation (10 1, and will therefore set the normalisation of 



our baryonic models. 



It is to be noted that equation ( 10 1 is obtained un- 



der the assumption that all the sources lie at the same 
distance Ds- This is no longer valid for microlensing 
observations towards the galactic centre, as self-lensing 
(caused by the fact that a "source" in a lensing event is 
also acting as lens for another event) plays an important 
role. In this case one has to integrate not only on the dis- 
tance of the lenses but also on the distance of the sources, 
since the number density of both sources and lenses can 
vary substantially [15! : 



(r)(£,6) = 



J^^dDsT{e,b,Ds)dns/dDs 
J^°° dDs dus/dDs 



(11) 



where dug/dDs ex ps{Ds)D'^^^^'' is the distance distri- 
bution of detectable sources, Ps is the mass density of 
sources and Too is the maximum distance at which sources 
can be found that we fix to Too = 20 kpc (our results do 
not depend much on the specific choice of this parame- 
ter). Notice that the expression for drig/dDg, introduced 
in [T5], includes both a volume effect in the D^ factor 
since more sources lie at greater distances, and a lumi- 
nosity effect in the factor Dl^" given that only sources 
above the threshold luminosity can be detected (assum- 
ing that the distribution of sources above luminosity L 
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FIG. 1: Measured and expected microlensing optical depths 
for different galactic regions. The measured optical depths 
along with the corresponding uncertainties are shown by the 
boxes, while the expected values from our baryonic models are 
shown by the blue dots and error bars. All five models defined 
in Section PI have been scaled to agree with the MACHO |23| 
measurement towards {£,b) = (1.5°, —2.68°) represented by 
the leftmost box. Using this normalisation we show the ex- 
pected optical depth for different fields of view observed by 
several campaigns: OGLE-II GB [H], EROS GCl-4 pS] and 
EROS eUus, 7Nor, 7Sct, ^Sct [26]. In this plot all uncertain- 
ties are la and Rq — 8 kpc. 



follows L^=). The index (3^ depends, of course, on the 
type of source and it is usually taken to be /3s < —1 for 
main sequence stars and /3s = for red giant clump stars 
[48l l63] . In the following Section we explain how we use 
Ps and pi to recover the expected optical depth from each 
of the models defined in Section Inj while comparing it 
with observations. 



B. Observations 

Early measurements of the optical depth towards the 
galactic bulge dH |17| |33] showed values significantly 
higher than what expected from initial estimates 14, l5]. 
This motivated the study in |32j . as previously discussed, 
in the then-justified belief that the higher density of 
lenses for microlensing would leave less room for the gas 
and Dark Matter components in our Galaxy. However, 
since those early measurements, huge progress has been 
done in microlensing observations and modelling of the 
bulge (see Section In] for this latter point): in fact, the 
MACHO, EROS and OGLE collaborations have recently 
performed measurements of the optical depth towards 
the galactic centre. The present values are less constrain- 
ing than the early results, thus making it compelling to 
perform a more accurate analysis of the different mass 
components of our Galaxy. 

As a default, we start by adopting the 2005 measure- 
ment of the MACHO collaboration towards the galactic 
bulge [531. Out of the 94 observed fields, they identify 9 



at very small galactic latitude and longitude as the Cen- 
tral Galactic Region (CGR), and find the average value 



2.17to:38 X 10"^ for {£,b) 



(1.50°, -2.68°) 



(12) 

We scale our bulge models to meet this optical depth 
measurement as done in |41) : assigned all model param- 
eters but po_b for models 1-4 or /o,6 for model 5, we set 
the latter by requiring that the optical depth obtained 
using equation (11) matches the value in equation (12 1. 



We checked that our optical depth results for models 1 
and 2 are in agreement with those obtained in Figure 1 of 
Ref. 41J. Since the source stars monitored in 23 are red 
clump giants in the galactic bulge, we fa. /Sg — in equa- 
tion (11) while using ps = pj, and pi = Pb + Pd- Notice 



that this procedure gives an la range for po^b or /o,t, cor- 
responding to the experimental uncertainty in equation 

We repeat the procedure for each model discussed in 
Section In] thus fixing the galactic baryonic component. 
It is now possible to cross-check our models against cur- 
rent observations from other experimental groups. We 
do so by computing the optical depth through equation 
( 11 ) making use of the just-defined bulge normalisations. 



In particular, we consider the following measurements: 



2.55 



+0.57 
-0.46 



10"^ for 



• OGLE-II GB [21]: (r) 

(i, b) = (1.16°, -2.75°). In this case we use (3^ = 0, 
Ps = Pb and pi = pb + Pd in equation ( 11 ) since the 



source stars are red clump giants in the galactic 
bulge. 

EROS GCl-4 [55]: (r) = 2.90 ± 1.30 x 10"^ for 
{£,b) = (-4.5°, 2.40°); (t) = 2.32 ±1.73 x IQ-^ for 
{£, b) = (-1.5°, 2.42°); (t) = 2.20 ± 1.56 x 10"^ for 
{l,b) = (1.5°, 2.22°); (r) ^ 1.65 ± 0.83 x 10"^ for 
{£,b) = (4.5°, 2.53°). In this case we use (3^ = 0, 
Ps = Pb and pi = pb + Pd in equation ( 11 ) since the 



source stars are red clump giants in the galactic 
bulge. 



t) = 0.67l[J-^^ X 10-6 for 
■ 0.49l°:?ix 



EROS spiral arms 

i£,b) = (306.56°, -1.46°) 6'Mus; (r 

10-<^ for {£,b) = (331.09°, -2.42°) 7 Nor; (t) = 

0.72+°-^^ X 10-6 for {£,b) = (18.51°, -2.09°) 7 Set; 

-1+0.23 



T) = U.5UA5:i5 X iO-6 for {£,b) = (26.60°, -2.15°) 
/? Set. In this case we use pi = pb + pd in equation 
(111 and replace dUs/dDg by a gaussian of mean 
7 kpc and standard deviation 1 kpc in good agree- 
ment with the source distance distribution shown 
in Figure 11 of Ref. [15]. 

The results of this cross-check are presented in Figure 
[1] where the boxes represent the experimental measure- 
ments and corresponding Icr uncertainties, while the blue 
dots and error bars refer to the theoretical expectations 
of the five baryonic models once normalised to the MA- 
CHO CGR result in equation (12). It is evident that we 



obtain a very good agreement for all models, since the 
expected ranges fall nicely within less than la off the 



measurements, the only exceptions being models 1 and 2 
for the spiral arm /3Sct (but in this case the discrepancy 
is just slightly above la and so we shall not consider it 
worrying). This is in fact a reassuring result. 

Notice that we reproduce the experimental results not 
only in the innermost galactic fields reached by OGLE-II 
[H] and EROS [2S], but also beyond the limits of the 
central galactic regions where EROS collected microlens- 
ing events from four spiral arms across the sky [21]. It 
is hence fair to state that the five models defined in Sec- 
tion ]lT] and normalised to 2005 MACHO optical depth 
in equation ( 12 ) are in good agreement with present mi- 



crolensing data. Moreover, we have further checked that 
the corresponding bulge masses are well within the values 
found in literature Mf, ~ 1 - 2 x 10^" Mq ;35',1I1[M1I1S]. 
We shall use these five models in the remainder of the 
present work. 

Several comments are in order here. First, let us stress 
that all microlensing data used up to now (including 
the 2005 MACHO result [13], which is the one adopted 
throughout our analysis) are sensitive to Einstein times 
3 ^ ^fi/days < 700, or in lens masses 1.4 x 10~^ < 
Mi/M© < 79 (in this simple estimate we have assumed 
the typical values vt = 100 km/s, Dg = 10 kpc and 
Di = Q.hDs). This means that compact, baryonic ob- 
jects of virtually all masses (including all stars and the 
most massive planets) contribute to the observed optical 
depth to which we are normalising our bulge models. 

In particular, less massive planets (Afj/M© < 1.4 x 
10-^) could also contribute significantly to the galactic 
mass budget, and one might wonder whether the popula- 
tion of unbound Jupiter-mass objects recently discovered 
through microlensing |66j could lead to a systematic un- 
derestimate of the baryonic model normalisation. The 
presence of these objects has been invoked to explain 
MOA and OGLE microlensing events with ts < 2 days, 
and their abundance has been estimated to be approx- 
imately twice that of main-sequence stars |66| . Their 
contribution to the mass density in the bulge/bar and 
disk is 



/pL — 



/ dM M{dN/dM)pL 
JdMM{dN/dM)st.r. 



(13) 



where {dN/dM)pi^ and {dN/dM)stars are the mass func- 
tions of the Jupiter-mass population and the stars, re- 
spectively. Taking a Dirac delta for {dN/dM)pi^ and 
model 1 of Table S3 of the supplementary information 
of Ref. [66], we find /pl ~ 0.1%. We shall therefore 
neglect this newly discovered population in our calcu- 
lations, since its contribution is much smaller than the 
uncertainties in the mass model. 



IV. ROTATION CURVE 

This Section is dedicated to outlining the formalism 
and observational status regarding one key dynamical ob- 
servable, the galactic rotation curve. 
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FIG. 2: Tlie brealc-down of the rotation curve Vc{r) for different Millcy Way mass models. Tlie lower solid, dotted and dashed 
lines correspond to the individual contributions of bulge/bar, disk and halo, respectively. In all frames the adopted Dark Matter 
halo follows a spherical NFW profile with r^ — 20 kpc, a = 1 and po = 0.4 GeV/cm . The rotation curve for model 2 has been 
omitted since it is very similar to the one featured in model 1. The upper blue solid line with error bars indicates the expected 
total rotation curve, while the red (dark) boxes show the compilation of data in [30] rescaled to i?o = 8 kpc and vq — 230 km/s. 
The leftmost shaded area shows the cut r > 2.5 kpc used throughout the analysis. The red long-dashed line in the bottom 
right frame shows the adiabatic compression of an initial NFW profile with r^ = 20 kpc, a — 1 and po = 0.4 GeV/cm and 
assuming model 5 for the baryonic mass distribution. In these plots all uncertainties are la. 



A. Theoretical framework 

Regardless of their nature, all forms of matter con- 
tribute to the rotation curve of our Galaxy, unlike the 
case of microlensing, where only compact bodies along 
the line of sight play a role. This is basically what allows 
us to extract information about the DM and gas com- 
ponents. In full generality, the circular velocity Vc at a 
given galactocentric distance r reads 



vUr) 



E<- 



w = E^7r("'^ = ^/2'^) 



(14) 



where i runs over all the mass components (bulge/bar, 
disk and Dark Matter halo) , 0i is the gravitational poten- 
tial associated to component i and (r, 9, ip) are galactic 
spherical coordinates {9 — 7r/2 defines the plane of the 
galactic disk). For the spherical Dark Matter halos in 
equations (|6| and ^ one recovers the well-known ex- 
pression vljjj^,j{r) = GMdm{< r)/r, where Mdm{< r) 
is the DM mass enclosed in a sphere of radius r around 
the galactic centre. The case of an oblate halo referred to 



at the end of Section |TT] can also be treated analytically 
with a slightly more complicated expression (cf. equa- 
tion (2-91) in [67). As for the baryonic contribution, 
let us notice that all components of model 5 are speci- 
fied through their gravitational potential, which renders 
trivial the computation of the corresponding circular ve- 
locity with equation (14 1. Finally, since the mass dis- 



tributions of models 1-4 are rather complicated and tri- 
axial in general, there is no simple expression for Vc] in 
this case we compute the gravitational potential due to 
an arbitrary mass distribution by expressing the solution 
of the Poisson equation as a series of multipoles up to 
order Imax — 2 (cf. equation (2-122) in [ST]) and then 
apply equation ( [l4| ). All theoretical expectations for Vc 
presented in the following have been averaged over ip in 
order to ease comparison with experimental data (which 
refer to different positions across the galactic plane). 

With this formalism at hand we can finally compute 
the rotation curve associated to the different models spec- 
ified in Section [n] and whose bulge/bar components were 
appropriately normalised to microlensing data (see Sec- 
tion III and Figure IT] for details) . Figure [2] displays the 



rotation curve break-down for our baryonic models (mod- 
els 1 and 2 yield very similar rotation curves despite the 
different bulge shapes, so we omit the latter model for 
plotting purposes) and selecting a spherical NFW profile 
with r, = 20 kpc, a = 1 and po = 0.4 GeV/cm^ The 
bulge/bar, disk and halo contributions are represented by 
the lower solid, dotted and dashed curves, respectively, 
in the range r = — 10 kpc with bins of 0.5 kpc. Notice 
that for model 5 the disk contribution in dotted contains 
both the stellar and gas components. The upper blue 
solid curve with error bars in every frame of Figure [2] de- 
notes the total rotation curve predicted by each model. 
The error bars show the propagated uncertainty due to 
the bulge normalisation scaled to MACHO microlensing 
optical depth in equation (12). As we shall see, the com- 



bination of microlensing data and the rotation curve sets 
an upper limit on the distribution of Dark Matter across 
the Galaxy. Before deriving these constraints, we turn to 
the status of present data on the galactic rotation curve, 
represented by the red (dark) boxes in Figure [2J 



B. Observations 

Different methods are available to infer the galactic 
rotation curve at galactocentric distances r < 10 kpc, 
which is the range we are interested in for the scope 
of this paper. One of the most used techniques relies 
on the observation of gas clouds moving in the galactic 
plane. For each line of sight {£, b) the extreme veloc- 
ity vt - the so-called "terminal velocity" - is measured 
and converted to the circular velocity at a given radius 
under the assumption of strictly circular gas orbits |68j : 
VciRosinl) — vt{t) + Wo sin ^, where vq = Vc{Ro) is the 
local circular velocity. The 21 cm line is widely used in 
literature to determine the terminal velocity of gas clouds 
in the inner regions of the Milky Way. Other kinematical 
probes for the outer regions of our Galaxy are C-stars, ob- 
served through near-infrared photometry (e.g. |69j). and 
the CO emission line of HII regions (e.g. [70]). 

Generically, the observed velocities of either gas clouds 
or stars with respect to the Earth must be transformed 
to a circular velocity through specific assumptions on the 
local reference frame. The different values adopted by 
the several observational groups over the years make it 
difficult to simply take their final values altogether. In 
Ref [30j the authors construct a rotation curve of the 
Galaxy starting from the very observational data, and 
unify the inference of rotation curves for a single set 
of local galactocentric radius and local circular velocity 
{Ro,vq) — (8.0 kpc, 200 km/s). Here we use that compi- 
lation of data, which includes measurements from several 
observational groups and techniques. 

The values of {Ro,vo) adopted in the literature vary 
significantly, and since the constraints we shall obtain 
in Section |V| depend quite sensitively on the choice of 
the Earth reference frame, we will take into account the 
error on Rq and vq throughout the analysis in order to 



"bracket" the uncertainty on the constraints. In addi- 
tion to the range Rq = 8.0 ± 0.5 kpc mentioned in Sec- 
tion |TT] we consider local circular velocities in the in- 
terval Wo — 230 ± 30 km/s. This choice reasonably en- 
compasses a broad set of measurements [71 - :75j . Recent 
analyses (e.g. [73]) in particular seem to indicate local 
velocities well in excess of the 1985 lAU recommended 
value Vq = 220 km/s [37 . In order to convert the data in 
[5U] to a given set {Rq,v'q), the following transformation 
laws are applied: 



R' = R 



Rq 



Vt 



R^ 

Rn 



K - vo) 



(15) 



with the primed quantities being the ones obtained with 
our new choice of {R'q, Vq), and the plain ones being those 
as from 30^, with (i?o,wo) = (8.0 kpc, 200 km/s). These 
transformation laws are only strictly valid for the mea- 
surement of terminal velocities. However, since our con- 
straints will come essentially from radii within the Solar 
circle, where terminal velocity data points dominate the 
compilation in Ref. [30j . this procedure is reasonable. 

Given that we are using a compilation of data obtained 
with different techniques, the treatment of uncertainties 
on the derived quantity Vc{r) is a delicate matter. The 
following strategy is adopted in this work. After scal- 
ing the data from [30] to a given (i?o,wo) according to 
equations (15), a binning of 0.5 kpc is applied in the 
galactocentric distance range r = — 10 kpc. Since sev- 
eral measurements lie in each r-bin, the central value on 
Vc is taken to be the mean of all the central values in the 
bin, while the corresponding la systematic uncertainty 
is the standard deviation of the central values. We es- 
timate the 1(7 statistical uncertainty as the mean of the 
la uncertainties of all the measurements in each bin, and 
add it to the systematic error. It is this total uncertainty 
that is represented by the red (dark) boxes in Figure [2] 
(for i?o = 8 kpc, Wo = 230 km/s) and that we use to 
gauge our knowledge on the Milky Way rotation curve in 
the following Section. 

We restrict our analysis to galactic radii r > 2.5 kpc, 
corresponding to longitudes |sin€| > 0.30 — 0.35 in line 
with what is usually done in the literature |431 165) . This 
cut - marked in the plots of Figure [2]- is imposed by the 
severe triaxiality of the gravitational potential in the in- 
ner galactic regions caused by the presence of the galactic 
bar, invalidating the assumption of circular orbits. 



V. RESULTS 

We now turn to the discussion of our results. From 
Figure [2] one can see how the modelling of the baryonic 
component in our Galaxy affects the rotation curve. First 
of all, let us recall that the expected circular velocity de- 
picted in that Figure is obtained by tuning our bulge 
models to microlensing data, independently of the data 
on the rotation curve itself. It is thus remarkable that 
in all frames the expectation roughly matches the data - 
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FIG. 3: Constraints on the Dark Matter distribution parameters po and a provided by current data on the Milky Way rotation 
curve for a generahsed NFW (left) and an Einasto (right) profile. The thick solid, thick dashed, thick dotted, thin solid and 
thick dot-dashed lines are the 2a constraints in the case where the galactic baryonic component follows model 1, 2, 3, 4 and 
5, respectively, all scaled to match the microlensing optical depth in equation (12 1. The green (light) shadowed area delimits 



the uncertainty on the constraint given the present-day baryonic models, while the red (dark) shadowed region indicates the 
excluded parameters. In the left frame, the upper thick solid line labelled "g = 0.7" refers to the constraints for an oblate NFW 
profile with axis ratio g = 0.7 in the case of taking model 1 for the baryonic component. The shadowed rectangle encompasses 
the ranges of profile slopes found in numerical simulations and the values of po found in the recent literature (see Section Il| , 
while the red filled circle in the left frame marks the parameter set (po = 0.4 GeV/cm , a — 1.0) used to produce Figure l2 In 
both frames we have fixed r^ = 20 kpc, i?o = 8.0 kpc and v^ = 230 km/s. 



this argues in favour of the accurateness of present Milky 
Way mass models in explaining both microlensing and 
dynamical observables. Models 1-4, in particular, suc- 
ceed in reproducing a flat rotation curve for r > 2.5 kpc 
but appear slightly below the data at 6 kpc ^ r < 10 
kpc, while model 5 does feature the exact same pattern 
indicated by the data. Again, we shall not regard the 
inner 2.5 kpc to set our constraints. 

Notice as well that the coarse r— binning may be mis- 
leading at small scales; model 5, for instance, produces a 
steep bump in Wc at r < 0.5 kpc (see e.g. Ref. [55]) which 
is not seen in Figure [2] simply because of the chosen bins. 



A. Conservative constraints on (poiQ) 

Figure [2] assumes an NFW profile with fiducial values 
Po = 0.4 GeV/cm^ a = 1.0 (and r^ = 20 kpc). We are 
now interested in studying the constraints that can be set 
in the DM parameter space (poj Q^) by using the data on 
the Milky Way rotation curve. A given set (poi «) will be 
considered to be excluded at 2cr if the la lower end of the 
expected Vc (curve with error bars in Figure pi) exceeds 
the 2(7 upper end of the measured Vc (boxes in Figure KJ) 
in at least one radial bin in the range 2.5 kpc < r < 10 
kpc. This is a rather conservative procedure indeed ~ 



note, for instance, that for model 4 the fiducial NFW 
setup (see bottom left frame of Figure [2]) is not excluded 
at 2(7, even though the predicted rotation curve is far 
from explaining the data. 

The constraints obtained with these definitions are 
shown in Figure |3] for the different models 1-5 and the 
generahsed NFW (left) and Einasto (right) DM pro- 
files. Red (dark) shading signals the excluded regions and 
green (light) shading spans the uncertainty on the exact 
constraint given the present-day baryonic models. These 
plots assume the fiducial setup r^ = 20 kpc, Rq — 8.0 
kpc, Wo = 230 km/s. It is worth mentioning that the con- 
straints presented in Figure [3] are generally dominated by 
two radial bins: the one centred at r = 2.75 kpc for small 
local DM densities and the one at r = 7.75 kpc for larger 
densities. 

The conservative constraints shown in Figure [3] 
are found to be consistent with the range of inner 
slopes/shape parameters a found in numerical simula- 
tions and local DM densities po hinted by dynamical ob- 
servables (see Section In] and shadowed rectangles in Fig- 
ure [3]). Hence, contrary to the findings of Ref. [35], we 
conclude that microlensing observations and data on the 
rotation curve of the Milky Way are perfectly compati- 
ble with the findings of iV— body simulations in ACDM 
cosmologies. Nonetheless, interesting constraints may be 



placed on the slope of the DM profile in our Galaxy: for 
the fiducial density po = 0.4 GeV/cm , Figured excludes 
a > 1.5 {a < 0.06) for the generalised NFW (Einasto) 
profile. Notice that in the NFW (Einasto) case one can 
set an upper (a lower) limit on a - this is simply be- 
cause for r <§:_ Ts, d\ogpDA[/d\ogr — —a for NFW and 
d\og pdm / dlogr = —2{r/rg)°' for Einasto. Therefore, 
unlike in the NFW case, a larger a for the Einasto pro- 
file corresponds to a faster roll of the slope to as r — > 
and thus a less steep profile. 

Deviations from spherical symmetry can in principle 
affect the constraints. By modelling an NFW oblate halo 
as detailed in Section lllj we obtain the exclusion curve 
labelled "q=0.7" in Figure [S] (left) for the case of model 
1: at first glance this constraint appears weaker than 
the spherical one, but it should be noted that an oblate 
profile corresponds to a higher po (about 20% higher ac- 
cording to |54j ) . Departures from spherical symmetry are 
therefore not able to weaken significantly our constraints. 

Up to now we have fixed the scale radius, galactocen- 
tric distance and local circular velocity to the respective 
fiducial values, r^ = 20 kpc, Rq = 8.0 kpc, vq — 230 
km/s. These astrophysical parameters, whose uncertain- 
ties are sizeable, affect in distinct ways our calculations. 
The scale radius r^, for instance, sets the concentration 
of the DM profile; the smaller rg the larger the DM con- 
tribution to the rotation curve. On the other hand, a 
smaller Rq shrinks the bulge and the disk leading to an 
increase in the bulge central density to produce the same 
optical depth; however, a smaller Rq also leads to a less 
constraining Vc data set so that overall the larger Rq 
the more aggressive our DM constraints. The local cir- 
cular velocity Vq, instead, sets essentially the plateau of 
the rotation curve and thus tighter constraints result for 
smaller vq. 

In view of these considerations and usin g th e ranges for 

■ ' -15 



2.0 



Rq and vq outlined in Sections II and IV (r^ = 20 



kpc, Ro — 8.0 ± 0.5 kpc, vq = 230^30 km/s), we define 
three astrophysical setups: (i) conservative, with rg = 35 
kpc, Rq = 7.5 kpc, Vq = 260 km/s; (ii) mean, with 
Ts = 20 kpc, Rq = 8.0 kpc, vq = 230 km/s; and (in) 
aggressive, with r^ = 10 kpc, Rq — 8.5 kpc, vq = 200 
km/s. The mean configuration was used in Figures ^ 
and [3j Figure |4] shows the effect of adopting the con- 
servative or aggressive setups on the derived DM upper 
limits for the generalised NFW profile. For simplicity 
we only show the upper limits encompassed by all mod- 
els, instead of individual constraints. We see from this 
Figure that, for reasonable local DM densities, an NFW 
profile in line with the findings of numerical simulations 
can only be (barely) excluded at the expenses of push- 
ing some astrophysical parameters to somewhat extreme 
values (in particular vq — 200 km/s). We are thus led to 
the conclusion that the results of Ref. p2] do not hold, 
given the available microlensing and dynamical data and 
our present knowledge on astrophysical parameters such 
as Ts, Rq or Vo- 
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FIG. 4: The bracketing of the 2a upper limits on the Dark 
Matter distribution parameters po and a for the generalised 
NFW profile and three astrophysical setups: conservative 
(dashed; r^ = 35 kpc, 7?o = 7.5 kpc, vo = 260 km/s), mean 
(solid; Ts = 20 kpc, Ro = 8.0 kpc, iiq = 230 km/s) and aggres- 
sive (dotted; r^ ~ 10 kpc, Ro = 8.5 kpc, vo — 200 km/s). The 
two lines for each setup encompass the upper limits set using 
the baryonic models 1-5. In particular, the mean shadowed 
area as well as the shadowed rectangle are the same as in the 
left frame of Figure [3] 



B. Determination of (po,q) 

Models 1-4 include no gas component, which in prin- 
ciple makes the corresponding exclusion curves in Figure 
|3] conservative. Model 5, instead, includes a gas disk and 
therefore we can go one step further and ask which DM 
parameters (po,a) provide the best fit to the measured 
rotation curve. 

For the fitting procedure we bin the Vc data as dis- 
cussed in Section IV A (0.5 kpc bins) and consider the 
radial range r — 2.5 — 10 kpc. The microlensing un- 
certainties were propagated through the baryonic model 
adopted, and the sum of statistical and systematic was 
used for the rotation curve data according to the pro- 
cedure described in Section HVBI We have assumed that 
different bins are uncorrelated: for the rotation curve ob- 
servations, different datasets dominate the measure (and 
the uncertainty) at different bins. While this reduces the 
correlation between different bins, it is also to be taken 
into account that the adoption of a systematic in addition 
to statistical is in line with our conservative approach. 

The results of this fitting procedure are shown by the 
contours (corresponding to Ax^ — 2.30, 6.18 for a two- 
parameters fit) in Figure p] for DM parameters inside 
the contours, model 5 manages to explain both the mi- 
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FIG. 5: Constraints on the Dark Matter distribution parameters po and a for a generalised NFW (left) and an Einasto (right) 
profile using the baryonic model 5. The thick dot-dashed curve is the 2ct constraint already shown in Figure |3J while the 
contours show the parameter space producing a good fit to the rotation curve {Ax^ = 2.30, 6.18) with the best-fit configuration 
indicated by the cross. The shadowed rectangle encompasses the ranges of profile slopes found in numerical simulations and the 
values of po found in the recent literature (see Secti on [HI ), while the red filled circle in the left frame marks the parameter set 
(po — 0.4 GeV/cm , a — 1.0) used to produce Figurel2rThe empty up-triangle, circle and down-triangle in the left frame show 
the local density and shape of the DM profile upon adiabatic contraction of the initial profile indicated by the corresponding 
filled symbols. The adiabatic contraction was applied using model 5 to fix the baryonic distribution Mi,{< r), that entails 
fb ~ 5.2%, 4.0%, 3.0% for the up-triangle, circle and down-triangle, respectively. In both frames we have fixed Vs = 20 kpc, 
Ro = 8.0 kpc and vg = 230 km/s. 



crolensing optical depth towards the bulge and the rota- 
tion curve of our Galaxy. It is interesting (and reassur- 
ing) that the contours fall nicely on top of the expected 
ranges indicated by the shadowed rectangles. In partic- 
ular, we find 



Po 
Po 



0.20 - 0.56 GeV/cm^ (NFW) 



(16) 



0.22 - 0.55 GeV/cm^ (Einasto) (17) 



locities, which dominate the compilation of Ref. [30] for 
r < Rq. In the range Ro < r < 10 kpc, the dominant 
data points on Vc{r) are due to the measurement of ve- 
locities in CO and HII regions. In that case the scaling 
for a given set {R[),v'q) is different from (151 and reads 
v' ~ V + {v'q — vq) a Rq = Rq = 8 kpc. This means 
that up to now we have been overestimating the scaled 
velocities in the range 8 kpc < r < 10 kpc by < 3% for 
Vq = 230 km/s. We have however checked that applying 



V' = V + (v'q — vo) for 8 kpc < r < 10 kpc and the set of 
transformations ( 15 ) for r < 8 kpc changes the la deter- 



at la, which is consistent with previous estimates ob- 
tained with different techniques [131 [S3] . 

We have checked the robustness of the result in equa- 
tion ([T6| in a number of different ways. First, the po = 0.17 - 0.50 GeV/cm'' in the case of NFW profiles 



mination of po by a small (but non-negligible) amount: 



role of the number and width of the radial bins in the 
fitting of the rotation curve was analysed. Increasing 
the size of the bins to 1 kpc leads to a measurement 
Po = 0.16 — 0.65 GeV/cm in the case of a generalised 
NFW profile, i.e. an increase in the uncertainty on po 
from 47% to 60% with respect to equation ( 16 1. A smaller 



binning would lead to a smaller uncertainty on po, but 
the average uncertainty on the galactocentric position of 
the rotation curve trackers in Ref. [3D] is ~ 0.5 kpc, which 
therefore provides a physical lower limit on the size of the 
bins. 

Second, as pointed out in Section [IV B[ the transfor- 
mations ( 15 ) are only strictly applicable for terminal ve- 



It is also interesting to see how the determination 
( 16 1 for NFW profiles changes for oblate halos. Us- 
ing the axis ratio q = 0.7, we obtain the la range 
Po = 0.26 — 0.68 GeV/cm . Using instead a more spher- 
ical halo with q = 0.9, our la local DM density determi- 
nation reads po = 0.22 — 0.59 GeV/cm , to compare with 
the recent measurement based on the vertical motion of 
stars in the solar neighbourhood [S51 [75]. Let us note 
that the up-shifting of the central value (with respect to 
the spherical case) is expected since the oblateness al- 
lows for the concentration of more DM along the galactic 
plane for an assigned halo mass. 
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We have also tested the adiabatic contraction model 
of Ref. [60] (outUned in Section [TT| and equation (Is])) in 
the case of the NFW profile. Using once again model 
5 for the baryonic component, we have contracted the 
initial profiles indicated in Figure [s] (left) by the filled 
up-triangle, circle and down-triangle with /b — Mh{< 
200 kpc)/Mtot(< 200 kpc) = 5.2%, 4.0%, 3.0%, respec- 
tively. The final DM profile turns out to be well fitted by 
a generalised NFW function with parameters marked by 
the empty symbols in the same Figure (the contracted 
profile corresponding to the filled circle is indicated by 
the red long-dashed line in the bottom right frame of 
Figure [2]). In particular, we find enhanced local DM den- 
sities and slopes a ~ 1.6 — 1.7, which are slightly above 
the value a = 1.5 found elsewhere [TT^ (see also refer- 
ences therein) but note that we are using the original 
adiabatic contraction model 6OJ and not one of its refine- 
ments [6TJ|62]. Although our analysis cannot rule out the 
presence of adiabatically compressed profiles since they 
depend on the initial total mass distribution and on the 
specific baryonic model adopted, it definitely allows us to 
claim that if the present-day DM profile is steeply rising 
towards the centre, then the local DM density must be 
small. For the specific case of a = 1.5 (1.7) we find an 
Icr range po - 0.25 - 0.34 (0.22 - 0.28) GeV/cni^. Some 
of the extreme models discussed in the literature, e.g. in 
the context of indirect DM searches |771 [75] , are therefore 
found to be ruled out by a combination of microlensing 
and dynamical observations. 



VI. CONCLUSIONS 

We have studied the constraints that microlensing and 
dynamical observations can set on the distribution of 
Dark Matter in the Galaxy, keeping into account all 
experimental uncertainties. Starting from state-of-the- 
art models for the galactic baryonic component, we have 
rescaled them to match the observed microlensing optical 
depth towards the galactic bulge, and compared the re- 
sulting rotation curve with the one inferred from terminal 
velocities of gas clouds and other kinematical probes. 

This allowed us to revisit the compatibility of different 
observational probes with the results that emerge from 
numerical simulations in ACDM cosmologies. We have 
followed two different approaches. In the first one, we 
have set conservative upper limits on the Dark Matter 
local density and profile shape towards the centre of the 
Galaxy, working with generalised NFW and Einasto pro- 
files. The fiducial parameters usually adopted in the lit- 
erature for both profiles have been found to be safely 
within the allowed regions set by our analysis, contrary 
to earlier claims of inconsistency between observations 



and cuspy Dark Matter profiles. 

In our second approach, we focussed on the only bary- 
onic model among those discussed here that also con- 
tains a description of the amount and distribution of 
gas, which is expected to provide a non- negligible con- 
tribution to the mass and therefore to the rotation curve 
of the Milky Way. For this specific model, we were able 
to calculate the values of the Dark Matter parameters 
that provide the best fit to the combination of microlens- 
ing and rotation curve data. The resulting Icr range for 
the local DM density, for NFW profiles, was found to 
be po — 0.20 — 0.56 GeV/cni , therefore consistent with 
estimates obtained in the literature with different tech- 
niques. 

Finally, we have studied the consequences of adiabatic 
compression, often invoked as a mechanism that could 
increase the amount of Dark Matter in the innermost re- 
gions of the Galaxy, and found that although our analysis 
is not able to discard this mechanism in general, it rules 
out combinations of local densities and profile slopes in- 
voked in the literature, e.g. in the context of indirect 
Dark Matter searches. For NFW profiles with a = 1.5 
(1.7) we constrain the local DM density to be within the 
Icr range pa ~ 0.25 - 0.34 (0.22 - 0.28) GeV/cm^. 

As numerical simulations that include the effect of 
baryons become ever more accurate and reliable, it will 
be interesting to compare the findings of these simula- 
tions with the wealth of observational data that are cur- 
rently available. In particular, these simulations should 
allow us to tie together the distribution of baryons and 
Dark Matter, and to provide a more precise prescription 
for the slope of the DM profile towards the galactic cen- 
tre. This slope is currently extrapolated from simulations 
without baryons, which is an unreliable procedure given 
that baryons are known to dominate the gravitational 
potential in the inner Galaxy. 
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